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Abstract 

An isotropic three-dimentional perfect lens based on cubic meshes of interconnected transmission 
lines and bulk loads is proposed. The lens is formed by a slab of a loaded mesh placed in between 
two similar unloaded meshes. The dispersion equations and the characteristic impedances of the 
eigenwaves in the meshes are derived analytically, with an emphasis on generality. This allows 
designing of transmission-line meshes with desired dispersion properties. The required backward- 
wave mode of operation in the lens is realized with simple inductive and capacitive loads. An 
analytical expression for the transmission through the lens is derived and the amplification of 
evanescent waves is demonstrated. Factors that influence enhancement of evanescent waves in the 
lens are studied and the corresponding design criteria are established. A possible realization of the 
structure is outlined. 
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I. INTRODUCTION 



In recent literature a lot of attention has been given to systems that are able to sense 
electromagnetic near fields (evanescent waves) and even to "amplify" them. The superlens 
proposed by Pendry 1 is one of such systems. His superlens is based on a Veselago medium 2 
slab. The real parts of the permittivity and the permeability of the Veselago slab are both 
negative at a certain wavelength. Thus, the eigenwaves in the slab are backward waves, 
i.e. the wave phase and group velocities are antiparallel. This provides negative refraction 
and focusing of waves in a planar slab, as was outlined by Veselago. 2 However, Pendry 
discovered that there was also a possibility to excite surface plasmon-polaritons on the slab 
surfaces and, due to that, amplify near fields 1 . The slab thickness can be of the order of 
the wavelength, so that the plasmon-polaritons excited at both sides of the slab are strongly 
coupled. Under certain conditions the plasmon-polariton excited at the back surface of the 
slab has a much stronger amplitude than that at the front surface. Such amplification of 
evanescent waves is the key principle in subwavelength imaging. 

Known experimental realizations of volumetric artificial materials with negative param- 
eters are highly anisotropic structures that utilize dense arrays of parallel thin conducting 
wires and variations of split-ring resonators 3 . Proposed isotropic arrangements use orthogo- 
nal sets of split rings 4 and also three-dimensional arrays of wires 5 . At the same time, there 
have been achievements in modeling of Veselago materials (and Pendry lens) with the help 
of LC-circuits or transmission-line (TL) based structures. 6-11 These networks do not rely on 
resonant response from particular inclusions, and the period of the mesh can be made very 
small as compared with the wavelength. These features allow realization of broadband and 
low-loss devices, which is extremely difficult if resonant inclusions are used. 

The transmission-line network approach has been successfully realized in one- and two- 
dimensional networks, but up to now there have been doubts if it is possible to design a 
three-dimensional (3D) circuit analogy of the Veselago medium. The difficulties arise from 
the fact that such a 3D network requires a common ground connector. Any realization 
of such a ground will effectively work as a dense mesh of interconnected conductors that 
blocks propagation of the electromagnetic waves practically the same way as a solid metal 
does (the structural period must be much less than the wavelength in order to realize an 
effectively uniform artificial material). In this paper we introduce isotropic three-dimensional 



2 



transmission-line networks that overcome this difficulty. 

In the TL-based networks that we study, the electromagnetic energy propagates through 
TL sections. The inside of every TL section is effectively screened from the inside of the 
other sections and from the outer space. This can be very naturally imagined with a 3D 
cubic-cell network of interconnected coaxial cable segments: The inner conductors of the 
segments are soldered at the network nodes; the same is done for the outer conductors. The 
whole system appears as a 3D pipe network where every pipe holds a central conductor and 
those conductors are crossing at the node points inside the pipes. A loaded TL network 
can be realized now by placing loading elements inside the "pipes". To couple the waves 
propagating inside the TL sections with the free-space waves one will have to apply a kind 
of antenna array with every antenna feeding a particular TL segment. 

When using transmission lines loaded with bulk elements we speak of waves in the mean- 
ing of discrete waves of voltages and currents defined at the loading positions. Let us note 
that in the TL sections as such the usual, forward waves propagate. Only because of the 
loading the discrete voltage and current waves appear as backward ones when appropriate 
loading impedances are used. 

While completing this manuscript, we learned about another possible design of a 3D 
transmission-line analogy of a backward- wave material described in Ref. 12. That design 
is based on Kron's formal representation of Maxwell's equations as an equivalent electric 
circuit. 13 In Ref. 12 only ID propagation was studied analytically and 3D properties were 
analyzed numerically. 

II. THREE-DIMENSIONAL TRANSMISSION-LINE NETWORKS 

The proposed structure of 3D super-resolution lens consists of two forward-wave (FW) 
regions and one backward-wave (BW) region. The 3D forward-wave networks can be realized 
with simple transmission lines and the 3D backward-wave network with inductively and 
capacitively loaded transmission lines. One unit cell of the BW network is shown in Fig. 1 
(the unit cell enclosed by the dotted line). In the 3D structure there are impedances Z/2 
and transmission lines also along the z-axis (not shown in Fig. 1). In view of potential 
generalizations, the loads are represented by series impedances Z/2 and shunt admittances 
F, although for our particular purpose to realize a backward-wave network, the loads are 
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simple capacitances and inductances. The unit cell of the FW network is the same as in 
Fig. 1 but without the series impedances Z/2 and shunt admittance Y. The equations that 
will be derived for these structures can be used in various implementations, but this paper 
will concentrate on the case when Z = 1/juC and Y = 1/jooL. 
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FIG. 1: Unit cell of a 3D backward- wave transmission line network (enclosed by the dotted line). 
The transmission lines and impedances along the z-axis are not shown. Transmission lines have 
the characteristic impedance Zq and length d/2 (d is the period of the structure). 



III. DISPERSION IN BACKWARD-WAVE AND FORWARD-WAVE NET- 
WORKS 

A. Dispersion equations 

First we will derive the dispersion relation for a simplified 3D BW network, i.e. we will not 
take into account the transmission lines. Such approximation is possible at low frequencies. 
Without the transmission line segments this derivation is quite simple and can be done by 
summing up all the currents that flow to the node (a, (3, 7) and equating this sum with the 
current flowing to the ground (through admittance Y), see Fig. 1. The result is 

+ u a + u a 

,13— 1,7 ,,3,7-1 - 6C/ a)/3i7 ) — U a ^ a Y. (1) 
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We look for a solution of the form U a $ a = U$e ji>r (r = r x xo + r y yo + r z z Q ), and if we use 
q x = k x d, q y = k y d and q z = k z d, (1) can be reduced to 
1 
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-My 



+ e 



+ e 



+jqx 
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6) = y, 



or 



cos(g :c ) + cos(g y ) + cos(q z 



ZY 



+ 3. 



(2) 



(3) 



If we now insert Z = 1/juC and Y = 1/juL we get the dispersion relation for the 
LC-loaded network: 

1 



cos(q x ) + cos(q y ) + cos(q z 



+ 3. 



(4) 



2oo 2 LC 

Next we want to take the transmission lines into account. The effect of the transmission 
lines can be derived by first evaluating a part of the three-dimensional network as the one 
shown in Fig. 2 and deriving the relation between the current that flows into a node and 
the voltages of adjacent nodes. 
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FIG. 2: Part of the 3D structure with transmission lines. 



If current l\ is flowing towards node (a, ft, 7) and the current that goes into node (a, ft, 7) 
from the left is I z , then using the ABCD-matrix for a transmission line we get: 



Ut = A t U 2 + BJ 2 , 



(5) 




(7) 



h = C t U 2 + DJ 2 , (6) 

where 

cos(k d/2) jZ Q sm(k Q d/2) 
^jZ ' 1 sm(kod/2) cos(k d/2) 
ko in (7) is the wavenumber of waves in the transmission lines. From (5) and (6) we can 
solve I\ and I 2 as functions of U\ and U 2 : 

U 2 {B % C t -A t D % ) + D t U x 



B t 



(8) 



* = (») 



(10) 
(11) 



Similarly for I 3 and I z we get: 

_ _ U 4 (B t C t -A t D t ) + D t U 3 

Us - A t U 4 
B t ' 

Next we can derive two equations for the current flowing through the series impedance Z 
and solve U<i from both of them: 

h = Z~\U 2 - U S ) = (12) 

^ U2 ~ A + z-i Bt ' (13) 

I, = Z-> { U 2 -U a ) = U ^-W + ™ (14) 

^ U 2 ^ • (15) 

If we let U2 in both equations be equal, we can solve C/ 3 as a function of Lq and U4: 

= UiZ^Bt - U,(B t C t - D t A t )(A t + Z^B,) 
3 (D t + Z-iB t )(A t + Z-iB t )-Z-*B? ■ [ ) 

In order to derive an equation for I z [the current that flows into node (a, (3, 7) from the 
direction of node (a, (3, 7 — 1)] as a function of Lq and C/ 4 , we insert (16) into (11) and get 

= Z-'BtU, - (B t C t - D t A t )(A t + Z- l B t )U 4 _ A 
2 [(A + Z-iB t )(A t + Z~iB t ) - Z-*B*]B t B t 4 - 1 j 

If we use Lq = U aAl _i and Lq = U a ^ m then I z = S B wU a ,p,j-i + K BW U a ^ tl . Because of 
the symmetry we can derive the dispersion relation exactly the same way as in (1) - (4), 
and for the case Z = 1/jtuC, Y = 1/juL the result is 

cosfe) + cos(g„) + cos(q z ) = 1 3^^, (18) 



where 



. = J"C 

,BW (A + juCB t ) (A t + juCB t ) + u?C*B\ ' 



(19) 



K -(B t C t -D t A t )(A t + juCB t ) _ A ( , 

BW [(A + juCB t )(A t + juCB t ) + u 2 C 2 B?]B t B t ' 1 Uj 
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To derive the dispersion relation for the forward-wave network, we can use the equations 
derived for the backward- wave network letting C — > oo and L — > oo. This way we get 
from (19) and (20) the following equations for Spw an d Kyw'- 

1 



'FW 



B t (A t + D t y 



(21) 



B t C t - AA 1 A t 
^ FW = " A + A B t ~ B t - (22) 



From (18) we get the dispersion relation: 

cos(q x ) + cos(g v ) + cos(g z ) = -3^-^. (23) 

B. Typical dispersion curves 

Dispersion curves for backward-wave and forward-wave networks can be plotted if the 
values of the transmission line parameters and L and C are fixed. Let us choose the param- 
eters of the TLs and the lumped components as: L — 10 nH, C = 5 pF, d = 0.012 m (the 
period of the network), Zo,tl,bw = 85 Ohm, Z 0i tl,fw = 85 Ohm (characteristic impedances 
of the TLs). See Figs. 3 and 4 for examples of dispersion curves when a ^-directed plane 
wave is considered (i.e. q x — q y — 0). k = ojy^/c, where c is the speed of light in 
vacuum. Notice that the BW network supports backward-waves only in the region where 
0.32 GHz< / <0.98 GHz. Above that frequency band and the following stopband, the BW 
network works as a normal FW network until the next stopband appears. 
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FIG. 3: Dispersion curve for a backward-wave network. C = 5 pF, L = 10 nH, d = 0.012 m, 
e r = 2.33, Z ,tl = 85 Ohm. 
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By tuning the capacitance C (or inductance L), the stopband between the BW and FW 
regions shown in Fig. 3 can be closed, see Fig. 5a, where C = Cq = 4.151 pF. As can be 
seen from Fig. 5b, by changing the value of C from this "balanced" case, the stopband is 
formed either by moving the edge of the FW region up (C < Co) or by moving the edge of 
the BW region down (C > C ). 
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FIG. 5: Dispersion curves for the backward-wave network. L = 10 nH, d = 0.012 m, e r = 2.33, 
Z ,tl = 85 Ohm. a) C = 4.151 pF. b) Solid line: C = 3 pF, dashed line: C = 6 pF. 

IV. THE CHARACTERISTIC IMPEDANCES OF BACKWARD-WAVE AND 
FORWARD- WAVE NETWORKS 

Next the characteristic impedance of the backward-wave network is derived. If we assume 
that the interface between the two networks is in the center of capacitor C (see Fig. 2, where 
Z = 1/juC), then we can define the characteristic impedance as 

u 2 + u 3 u 2 + u 3 



j o,bw 



(24) 



2/ 2 2/ 3 

First we have to express U 2 , U 3 and / 3 (or optionally I 2 ) as functions of U\ and f/ 4 . We can use 
equations (14) - (16). If we insert (16) into (15) and (14), we find U2, U% and /3 as functions 
of Ui and £7 4 . Therefore we can present U 2 , U 3 and / 3 simply as: U 2 = obw^i + &bw^4i 
U 3 = cbw^i + 4w^4 and / 3 = e BW £/i + /bw^4- 

Because U4 = e~ jQz Ui for a wave moving along the +z-direction, the characteristic 
impedance can be expressed as 

U 2 + U 3 a BW + b BW e~ jqz + c BW + d BW e~ jqz 



J 0,BW 



2h 



where 



a B w 



2e BW + 2/Bwe-i* 
D t +jLuCB t 



(A + juCB t )(A t + juCBt) + uo 2 C 2 B 



2 • 



B t C t -D t A t D t +juCB t 



{B t C t -D t A t ){A t + juCB t 



]u;CB t 



ju;CB t (A + juCB t )(A t + juCB t ) + u*C*B? ' 



(25) 

(26) 
(27) 
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CBW — /ri x ■ i .nr>\rA , .r^T3\ , , e?2 ' 



(A + jwCA) (A + jwCA) + " 2 C 2 A 2 

(AA~ D t A t )(A t +jujCB t 
(A + j'^CAXA + jcuCA) + cu 2 C 2 A 2 
A jwCA 



2 ' 



(29) 



6BW A (A + juCB t )(A t + juCBt) + u 2 C 2 B? ' (30) 



A A - AA A (A A - A A) (A t + jcuCA; , . 

/mv ' A A(A + ^CA)(A + ^CA)+^ 2 C 2 A 2 ' 1 j 

To derive the characteristic impedance of the forward-wave network, we can use the 
equations derived for the backward-wave network if we insert C — > oo in them. If this 
condition applies, we get from (25) - (31): 

„ U 2 + U 3 a FW + b FW e~ jqz + c FW + d FW e~ jqz 



J 0,FW 



2/ 3 2e FW + 2/Fwe-"' 



(32) 



where 



a FW — c FW — — — — pr, (33) 
A + u \ 

Ofw = Afw = , , n , (34) 

A 1 

eFW = A ATA ' (35) 
AA - A A A A A - AA 

/fw " A A A + A • (36) 

Zo,bw and Z 0i fw can be plotted from (25) and (32) as functions of the frequency if the 
transmission line parameters and L and C are fixed. Let us choose the parameters of the 
TLs and the lumped components as: C = 5 pF, L — 10 nH, d = 0.012 m, e r . = 2.33, 
Zo,tl,bw — 85 Ohm, Zo,tl,fw = 70 Ohm (characteristic impedances of the TLs). See 
Figs. 6 and 7 for examples of the characteristic impedances when a z-directed plane wave 
is considered (i.e. q x — q y — 0). The effect of changing C on the characteristic impedance 
can be seen in Fig. 6, and the effect of changing Zo,tl,fw on the characteristic impedance is 
shown in Fig. 7. Notice that for the BW network the characteristic impedance is continuous 
only in the "balanced" case (C = 4.151 pF here), because in the stopbands the real part of 
the impedance is zero. 
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/ [GHz] 

FIG. 6: Characteristic impedance as a function of frequency for the backward-wave network. 
L = 10 nH, d = 0.012 m, e r = 2.33, Z ,tl = 85 Ohm. Dashed line: C = 3 pF, dash-dotted line: 
C = 4.151 pF, solid line: C = 5 pF, dotted line: C = 7 pF. 
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FIG. 7: Characteristic impedance as a function of frequency for the forward- wave network, d = 
0.012 m, e r = 2.33. Dotted line: ^o,tl = 50 Ohm, solid line: Zo,tl = 70 Ohm, dash-dotted line: 
■^o,tl = 90 Ohm, dashed line: ^o,tl = 110 Ohm. 

V. MATCHING OF FORWARD- WAVE AND BACKWARD- WAVE NETWORKS 

We consider a perfect lens with axis parallel to the z-axis. To have perfect imaging, the 
lens should support all spatial harmonics (i.e. waves with all possible transverse wavenum- 
bers k t ) of the source field and for those values of k t , /c 2j fw an d & 2jB w should be equal in 
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magnitude but opposite in sign. 

From Fig. 8 we can conclude that the matching of |/c 2j fw| and |/c 2) bw| (which corresponds 
to a relative refraction index of —1) can be achieved only at one frequency (depending on the 
parameters of the forward- wave and backward- wave networks). In Fig. 8 this frequency is 
/ = 0.7279 GHz. At this frequency the dispersion curves of the forward- wave and backward- 
wave networks intersect. In the analytical form this means that 

as can be seen from (18) and (23). The dispersion curves in Fig. 8 are plotted so that q x 
and q y are zero (a plane wave moving along the z-axis). From Fig. 8 it is seen also that 
for / = 0.7279 GHz | 

Qz,bw\ — |<7z,fw| = 0.4869 when k t = 0. This means that the absolute 
value of the maximum wavenumber for propagating waves is approximately 40.6 m _1 . 
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FIG. 8: Dispersion curves for the forward- wave (solid line) and backward- wave (dashed line) net- 
works. C = 5 pF, L = 10 nH, d = 0.012 m, e r = 2.33, Z ,tl,bw = 85 Ohm, Z ,tl,fw = 70 Ohm. 
The transverse wavenumber k t = 0. 

In addition to matching the wavenumbers (refractive indices), to realize an ideal "perfect 
lens" the interfaces between the forward-wave and backward-wave networks should be also 
impedance-matched. If the two regions would not be matched, reflections from the interface 
would distort the field patterns both inside and outside the lens. As can be seen from Figs. 6 
and 7, the characteristic impedances of the backward-wave network and the forward-wave 
network are about 40 Ohms at the frequency where the wavenumbers are matched (and 
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when C = 5 pF and -Zo,tl,fw = 70 Ohm). Notice that the impedances of the transmission 
lines in the forward-wave network have been lowered from 85 Ohm to 70 Ohm to achieve 
impedance matching of the forward-wave and backward-wave networks. 

Next the effect of nonzero k t on the matching is considered (at the optimal frequency). 
The minimum and maximum values of k t can be found from k tmin = —n/d and k t ^ max = n/d. 
From (18) and (23) we can plot /c^.bw and A; 2j fw as functions of the transverse wavenumber 
(k t = Jk% + ky) if we fix the frequency. Now /c 2j bw and k Zt pw are surfaces with variables 
k x and k y . By comparing these surfaces, it was seen that they are practically the same for 
all possible values of k t . This happens only at the frequency / = 0.7279 GHz, where the 
dispersion curves of the forward- wave and backward- wave networks intersect. 

Because the characteristic impedances are functions of the ABCD-matrices and k z (and 
k z is a function of k t ), we can plot Z 0i fw and Zo,bw as functions of k t if we fix the frequency. 
Now -Zo.bw and Zo,fw are surfaces with variables k x and k y . By comparing these surfaces, 
it was seen that they are almost the same (less than one percent difference) for all possible 
values of k t at / = 0.7279 GHz. See Fig. 9 for a 2D cut of relative difference between such 
surfaces (now k y = and therefore k t = k x ). The ideal "perfect lens" situation is achieved 
with Z 0) tl,fw = 70.59 Ohm, as can be seen also from Fig. 9. 
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FIG. 9: Relative difference (Zo,bw — -^o,fw)/-^o,bw as a function of k t . C = 5 pF, L = 10 nH, 
d = 0.012 m, e r = 2.33, Zo,tl,bw = 85 Ohm. Solid line: ^o,tl,fw = 70 Ohm, dashed line: 
■^o,tl,fw = 70.59 Ohm. k x = k t . 

When the frequency deviates from the optimal value (for which (37) is true), the 
wavenumbers and characteristic impedances are no longer matched. The effect of this is 
distortion of the image seen in the image plane of the lens. 
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VI. TRANSMISSION OF THE SOURCE FIELD THROUGH THE LENS 



The transmission coefficient of the lens can be solved by considering the incident and 
reflected fields in the lens system. Let us assume that the incident field outside the lens has 
the unit amplitude, the reflected field outside the lens has amplitude R, the incident field 
inside the lens has amplitude Ae~^ kz ' BwZ , and the reflected field inside the lens has amplitude 
^ e +jfe,Bw2 ( z j s the distance from the front edge of the lens). From these values we can form 
the following equations (the length of the BW slab is / and the transmission coefficient of 
the BW slab is T Lens , see Fig. 10): 



1 + R = A + B, (38) 




FW 

FIG. 10: A 3D perfect lens with distances I, s\ and S2 shown. 
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The resulting equation for the transmission coefficient liens is: 

rp (i \ 4Zo,FW^O,BW 
J- Lensl^t 



(42) 



(Z , FW + Z 0i Bw) 2 e+^- Bwi - I Z„ 
The total transmission from the source plane to the image plane (see Fig. 10) is then (distance 
from source plane to lens is S\ and distance from lens to image plane is s^) 



^tot(^t) — ^~Lens(^t) e 



-j&z,Fw(si+S2) 



(43) 



The longitudinal wavenumber k z as a function of k t can be found from the dispersion 
relations. Let us choose k t = k x and k y = so we can plot curves instead of surfaces. T tot 
as a function of k t can now be plotted if the frequency is fixed. Let us choose the lengths of 
the lens system as the following: / = 0.12 m, si = 0.06 m, s 2 = 0.06 m. Now we can choose 
the frequency at which we want to calculate T tot . 

Let us study the transmission properties at the matching frequency / = 0.7279 GHz. 
From (42) and (43) we can plot the magnitude and phase of T tot as a function of k t , see 
Fig. 11, case 1. From Fig. 11 it is seen that the "lens" works quite well for the propagating 
modes (—40 m _1 < k t < 40 m" 1 ), see Fig. 12 for an example of phase correction in the 
image plane. 
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FIG. 11: Total transmission from the source plane to the image plane as a function of kt- C = 5 
pF, L = 10 nH, d = 0.012 m, e r = 2.33, Z ,tl,bw = 85 Ohm, Z 0j tl,fw = 70 Ohm. Case 1: 
Sl = S2 = 0.06 m, I = 0.12 m. Case 2: s x = s 2 = 0.024 m, I = 0.048 m. Solid line: |T to t| (case 1), 
dashed line: arg(Ttot) (case 1), dash-dotted line: |T to t| (case 2), dotted line: arg(T to t) (case 2). 
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FIG. 12: Phase of electric field (propagating wave), a) Source plane is at z = m. b) Back edge of 
the lens is at z = m and image plane is therefore at z = 0.06 m. C = 5 pF, L = 10 nH, d = 0.012 
m, s r = 2.33, Zo,tl,bw = 85 Ohm, ^o,tl,fw = 70 Ohm. s\ = S2 = 0.06 m and Z = 0.12 m in T tot . 

According to Fig. 11, for evanescent modes (k t < —41 m -1 ,^ > 41 m _1 ) the "lens" 
works only in a limited range of k t , where the absolute value of the transmission coefficient 
T tot is greater than zero. One can notice that for evanescent modes a mismatch in k z 
affects mostly the phase of T tot in the area of propagating waves and a mismatch in the 
characteristic impedances affects primarily the absolute value of T tot . To improve the effect 
of evanescent waves enhancement, the characteristic impedances should be matched better 
in the evanescent wave area of k t (i.e. k t < —41 m _1 , k t > 41 m _1 ). There are several ways 
to achieve a better matching of the characteristic impedances. 

First, there is of course a possibility to change the impedances of the transmission lines, 
but this is probably not practically realizable because it would require very accurate manu- 
facturing (even a very small deviation from the ideal impedance values destroys the effect of 
growing evanescent waves). The tuning of Z 0i tl,fw was t es t e d and using the exact impedance 
required (see Fig. 9), the transmission of evanescent waves was clearly improved. The reso- 
nance peaks in Fig. 11 were moved further away from the center and the absolute value of 
T tot was larger than or equal to unity approximately for —100 m _1 < k t < 100 m" 1 ). 

Second, there is a possibility to change the frequency and study if the impedance matching 
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can be made better that way (this also means that the matching of wavenumbers k z is made 
worse which can also destroy the effect of growing evanescent waves). This was tested and 
the best results were obtained using frequency / = 0.72905 GHz. The region of transmitted 
fc t 's was again increased, i.e. the resonance peaks in Fig. 11 were moved further away from 
the center and the absolute value of T tot was larger than or equal to unity approximately 
for -100 m" 1 < kt < 100 m" 1 ). 

The third way to enhance the growth of evanescent waves is to change the length of 
the "lens". From (42) it is seen that the growth of evanescent waves is destroyed by the 
term (Z^fw — Zo,Bw) 2 e -: ' fcz ' BwZ in the denominator. This term can be made smaller by 
decreasing the length of the "lens" /. See Fig. 11, case 2 (|T tot | is larger than or equal to 
unity approximately for —160 m _1 < k t < 160 m -1 ) and Fig. 13, where the distances equal 
I = 0.048 m, s\ = 0.024 m and s 2 = 0.024 m. From Fig. 13 one can conclude that there is a 
significant growth of evanescent waves in the lens. 

1 

|E| 
0.5 





-0.03 0.03 -0.03 0.03 

x [m] x [m] 

a) b) 

FIG. 13: Absolute value of electric field (plots normalized to the maximum value) at / = 0.7279 
GHz. a) Solid line: source field (z = m), dashed line: field in the back edge of the lens (z = 0.072 
m). b) Dotted line: field in the front edge of the lens (z = 0.024 m), dash-dotted line: field in 
the image plane (z = 0.096 m). The source at z = produces only evanescent waves. The source 
consists of 10 harmonics with kt = (27r)/(10- 0.012 m) -n, where n = —5. ..5. C = 5 pF, L = 10 nH, 
d = 0.012 m, e r = 2.33, Z ,tl,bw = 85 Ohm, Z ,tl,fw = 70 Ohm, s\ = s 2 = 0.024 m, I = 0.048 
m. 
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By using the shortened lens and at the same time tuning the frequency appropriately, it 
was seen that the transmission coefficient could be made practically ideal (i.e. |T tot | = 1 
and arg(Ttot) = for all possible values of k t ). Using the shortened lens (same values as 
in Fig. 13) and frequency / = 0.7292 GHz, the absolute values of evanescent fields were 
indeed almost the same in the image plane and in the source plane (less than one percent 
difference) . 

VII. SUGGESTIONS FOR A PRACTICALLY REALIZABLE STRUCTURE 
A. Proposed structure 

How to manufacture three-dimensional transmission line networks? The main problem 
is the ground plane, which should exist in all three dimensions. One solution would be to 
use coaxial transmission lines (regular in the forward-wave network and loaded with lumped 
L- and C-components in the backward- wave network) as shown in Fig. 14a. This structure 
is realizable, but we propose a simpler structure based on microstrip lines, as presented in 
Fig. 14b. The problem with microstrip lines is of course the design of intersections where the 
transmission lines from six directions meet. This problem can be overcome by using ground 
planes which have holes in them at these intersection points. This way the conducting strip 
can be taken through the substrate and thus connection of the vertical conducting strips 
becomes possible. 




a) b) 

FIG. 14: a) Unit cell of a 3D forward- wave coaxial network, b) Unit cell of a 3D forward- wave 
microstrip network. L and C can be easily implemented as lumped components for both types of 
backward- wave networks. 
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B. Simulation results 



The proposed structure (see Fig. 14b) has been simulated in Ansoft HFSS (Version 9.2.1). 
Due to complexity of the structure and the limited calculation power available, the three- 
dimensional structure was simulated only near the first lens interface of Fig. 10. The simu- 
lated model had 10 x 3 x 3 (x x y x z) unit cells in the forward-wave region and 10 x 3 x 3 
unit cells in the backward-wave region. The properties of the transmission lines and lumped 
elements were the same as in Fig. 8. The edges of the system were terminated with matched 
loads to prevent reflections (Rbw — 85 Ohm in the backward-wave region and Rfw — 70 
Ohm in the forward- wave region). Different types of source fields (plane waves with dif- 
ferent incidence angles and a point source) were tested and in all cases negative refraction 
was observed at the interface between the forward-wave and backward-wave networks at the 
expected frequency (/ = 0.7279 GHz). 

A two-dimensional cut of the proposed structure was simulated as a complete "lens" 
system. Again negative refraction was seen at both interfaces, and therefore also focusing 
of propagating waves was observed. See Fig. 15 for the plot of the phase of the electric field 
in the two-dimensional simulation. The source field is excited at the left edge of the system 
in Fig. 15. When the field magnitude \E\ is plotted and animated as a function of phase, 
it is clearly seen that the phase propagates to the right in the forward-wave regions and to 
the left in the backward-wave region. The energy of course propagates to the right in all 
regions. 
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FIG. 15: Plot of &rg(E z ) in a two-dimensional part of the proposed lens system. The source field 
is a plane wave with incidence angle 9 = 25°. The source field is excited at the left edge of the 
system. C = 5 pF, L = 10 nH, d = 0.012 m, e r = 2.33, Zo,tl,bw = 85 Ohm, ^o,tl,fw = 70 Ohm. 

VIII. CONCLUSIONS 

In this paper we have introduced and studied a three-dimensional transmission-line net- 
work which is a circuit analogy of the superlens proposed by Pendry. The structure is a 
3D-network of interconnected loaded transmission lines. Choosing appropriate loads we 
realize forward-wave (FW) and backward-wave (BW) regions in the network. The disper- 
sion equations and analytical expressions for the characteristic impedances for waves in FW 
and BW regions have been derived. A special attention has been given to the problem 
of impedance and refraction index matching of FW and BW regions. From the derived 
dispersion equations it has been seen that there exist such a frequency at which the cor- 
responding isofrequency surfaces for FW and BW regions coincide. Theoretically this can 
provide distortion-less focusing of the propagating modes if the wave impedances of FW 
and BW regions are also well matched. Impedance matching becomes even more important 
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when the evanescent modes are taken into account. In this paper we have shown that the 
wave impedances can be matched at least within 1% accuracy or better if the characteristic 
impedances of the transmission lines are properly tuned. However, from the practical point 
of view an accuracy better than 1% becomes hardly realizable. It has been shown that 
decreasing the thickness of the BW region reduces the negative effect of the impedance mis- 
match, while the amplification of the evanescent modes is preserved. We have also outlined 
a couple of prospective designs of the perfect lens discussed in this paper and numerically 
simulated their performance. 
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